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>• ■ Abstract 

<N. 

Using coherent-state representations of quantum mechanics (Bargmann, Husimi, 
and stellar representations), we describe analytically the phase-space structure of 
the general eigenstates corresponding to a 1-dimensional bilinear hyperbolic Hamil- 
tonian, H = pq or equivalently H = \{P 2 - Q 2 ). Their semi-classical behaviour is 
discussed for eigenvalues either near or away from the separatrix energy {H = 0}, 
especially in the phase-space vicinity of the saddle-point (q,p) = (0,0). 

2 ■ 1 Introduction 

Semiclassical theory, and the quantum-classical correspondence, are still incom- 
pletely understood at the level of long-time or invariant structures, especially when 
the classical dynamics shows exponential sensitivity to initial conditions (instability, 
or positive Lyapunov exponents). 

We have therefore selected a simple 1-d system, the dilation operator which 
quantizes the classical linear hyperbolic Hamiltonian H{q,p) = qp, and collected 
some detailed analytical properties of its eigenstates. Elementary as this quantum 
Hamiltonian may be, seemingly unreported expressions for its eigenstates are given 
here, precisely within the coherent-state (Bargmann or Husimi) formulations where 
the semiclassical behaviour (K — > asymptotics) of quantum states is best seen. 

Formulae for this "dilator" should be useful tools to probe quantum phenomena 
linked to unstable (here, hyperbolic) classical dynamics. Here are two actively 
studied examples involving this Hamiltonian: 

- in one dimension, classical hyperbolic dynamics takes place near a saddle- 
point (assuming it is generic, i.e., isolated and nondegenerate) ; then it is locally 
equivalent to H = pq in suitable canonical coordinates, i.e, the classical dilator is 
the local normal form for this class of problems. At the quantum level, 1-d saddle- 
points also challenge semiclassical analysis: like energy minima, they correspond to 
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critical energy values (i.e., the phase-space velocity vanishes); semiclassical analysis 
is fundamentally harder at critical than at regular energy values (where simple WKB 
theory works), but saddle-points are even harder to understand than minima and 
their study has expanded more recently |IJ |2|, |3| . 

- in higher dimensions, the search for quantum manifestations of unstable classi- 
cal motion forms one facet of the "quantum chaos" problem. For instance, the role 
and imprint of unstable periodic orbits upon quantum dynamics remain actively de- 
bated issues. It is known that they influence both the quantum bound state energy 
spectrum (through trace formulas) and wave functions (through scarring). More 
specifically, the (stable) orbits of elliptic type generate quantization formulas and 
quasimode constructions in a consistent way. In chaotic systems, by contrast, trace 
formulas diverge and scarring occurs rather unpredictably, so that periodic orbits 
(now unstable) have incompletely assessed quantum effects; nevertheless, an essen- 
tial role must still go to their linearized classical dynamics, which is of hyperbolic 
type hence generated by linear dilator (s). 

Fully chaotic behaviour, a combination of global instability plus ergodic recur- 
rence, cannot however be captured by integrable models; therefore, the place of 
the 1-d dilator eigenfunctions in quantum structures corresponding to fully chaotic 
dynamics remains to be further studied. At present, these eigenfunctions should 
primarily find use as microlocal models for general 1-d eigenstates near saddle- 
points (and for separatrix eigenstates in higher-dimensional integrable systems, by 
extension) . 



2 Coherent-state representations 
2.1 Bargmann representation 

The Bargmann representation is a particular coherent-state representation of 
quantum wave-functions fj], in terms of entire functions. Although it can be 
defined in any dimension, we will just use it for 1-dimensional problems; it then 
transforms Schrodinger wave-functions ijj(q) defined over the whole real line into 
entire functions ip(z) of a complex variable 

= W> = TTuTI / e!H(* 2 +*' 2 )+^') jtf) dq>. (1) 

Here \z) denotes a (Weyl) coherent state localized at the phase space point (q,p) 
where z = 2~ 1 ^ 2 (q — ip), and satisfying (z\z') = e zz '/ h ; i.e. these coherent states are 
not mutually orthogonal, and their normalization is not {z\z) = 1, allowing instead 
the bra vector (z\ to be a holomorphic function of its label z. On the other hand, 
a closure formula exists which makes the Bargmann transform invertible, 

H = / e"^ \z)(z\. (2) 

Jc nfi 

The Bargmann transformation maps ordinary square-integrable wave-functions 
into a Hilbert space of entire functions of order < 2. We will however mostly 
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deal with generalized wave-functions ip(q), which are not I? but only tempered 
distributions. They can then still be Bargmann-transformed by the integral formula 
(1), and into entire functions of order < 2 as before, now bounded as 

i 2 |2 

\^{z)\ < c(l + \z\ 2 ) N evr for some N. (3) 
This in turn constrains the distribution of their zeros ||] : the counting function 

n(r) = #{zeros z m of -4>(z) s.t. \z m \ < r} (4) 
(zeros will be always counted with their multiplicities) verifies 

hmsup < (5) 

It follows that a Bargmann function admits a canonical Hadamard representation 
as an (in)finite product over its zeros, 

*(*)=ePM*»«» n (i_JL)exp(f + i(f) 2 ) (6) 

where p(z) = fC^ia^z 1 + aiz + ao), with moreover |a2| < 1/2; the integer n(0) > 
is the multiplicity of z = as a zero of VK- 2 )- This decomposition shows that the 
wave-function is completely determined by the knowledge of all the Bargmann zeros 
(with their multiplicities) plus three coefficients a2,a\, and ao (which only fixes a 
constant factor). 

When the Hamiltonian is a polynomial (or more generally an analytic function) 
in the variables (q,p), it can be useful to express its quantum version as a pseudo- 
differential operator acting on the Bargmann function 0, using the rules 

(creation operator) — ► multiplication by z 

d 

(annihilation operator) — ► h—. (7) 

dz 

2.2 Husimi representation 

An alternative point of interest lies in certain semi-classical densities on phase space 
associated to wave- functions. In particular, the Wigner function is defined as 

Vfy(g,p) = (2tt^)- 1 f ij}{q - r/2)^(q + r/2) e ipr ' h dr (8) 

and the Husimi function Q as the convolution of the Wigner function with a phase- 
space Gaussian, 

H f (q,p) = (nh)- 1 f W f (q',p')e-l^ 2+ ^ 2 V n dq'dp' . (9) 



+ 1 ~ W 
1 = 7=- 

q + ip 
a = — p— 
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The Wigner representation has greater symmetry (invariance under all linear sym- 
plectic transformations), but Wigner functions show a much less local semiclassical 
behaviour: they tend to display huge nonphysical oscillations, which must be av- 
eraged out to reveal any interesting limiting effects. In the Husimi functions, the 
spurious oscillations get precisely damped so as to unravel the actual phase-space 
concentration of the semiclassical measures, but at the expense of reducing the 
invariance group. The Husimi function is equivalently given by 

n^z) = = mz)\ 2 e-*/» (10) 

\Z\Z) 

hence it constitutes the density of a positive measure on the phase space; for the 
scattering-like eigenfunctions to be studied here, this measure will not be normal- 
izable. (For a normalized state, it is a probability measure thanks to the closure 
formula @.) 

It is interesting to study how the Husimi measure of an eigenfunction behaves 
as h — * 0. For an energy away from the separatrix, a standard theorem states that 
this measure concentrates on the energy surface, with a Gaussian transversal profile 
|io| , |lll , [l^l - For energies close to a separatrix, careful analyses were performed in 
|l|, ^, |||. Our aim here is to select a simple tractable case, namely the eigenfunctions 
of the linear hyperbolic Hamiltonian [pi], and to carry further its description by 



means of the Bargmann representation, using eq. (10) to derive the Husimi density 
as a by-product. 



2.3 Stellar representation 

According to the factorized representation (|6|), 1-d quantum wavefunctions can be 
essentially parametrized in a phase-space geometry by the distribution of their 
Bargmann zeros which, by eq. dlCf), is also the pattern of zeros for the Husimi 
density itself; it thus constitutes a complementary viewpoint to the previous em- 
phasis put on the high-density behaviour of the Husimi function. We refer to this 
"reduced" description of a wavefunction by a discrete cloud of phase-space points as 
a stellar representation. It puts quantum mechanics in a new perspective |l 
but calls for a finer understanding of both dynamical and asymptotic properties of 
Bargmann zeros if new results are to be awaited through eq. (^). Consequently, our 
subsequent analysis of "toy" eigenfunctions will largely deal with explicit behaviours 
of their Bargmann zeros. 



3 Description of the framework 
3.1 The linear hyperbolic Hamiltonian 

The classical 1-dimensional Hamiltonian of linear dilation is H(q,p) = pq, which is 
also equivalent to the scattering Hamiltonian H = \{P 2 — Q 2 ) upon a symplectic 
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rotation of the coordinates by tt/4 according to 

p + Q p-Q , in 

A classical trajectory at any energy E 7^ is a hyperbola branch, 

q(t) = q e\ p(t) = p e~* with E = p q (12) 

whereas the E = set is a separatrix, made of a stable manifold {p = 0}, an 
unstable manifold {q = 0}, and the hyperbolic fixed point (0, 0). 

We study eigenfunctions of the operator obtained by Weyl quantization, namely 
H = j(q-^ + \)- This quantum Hamiltonian admits two independent stationary 
wave-functions for any real energy E: 

^| = K9(±q) e^f-D 10 ^! (13) 

where 9{q) is the Heaviside step function; K ^ is a complex constant (having no 
preferred value, since the solutions are not square-integrable) . Microlocally, each of 
these wavefunctions is supported by the lagrangian manifolds A E = {pq = E , ±q > 
0}, i.e. half of the E'-energy surface Q. 

In order to obtain more semi-classical information, we will use the Bargmann 
representation. For instance, 

(z\^ E ) = ^(z) = J" eH-hW*^****) _L dq > (14 ) 

and we have of course 



^ E {z) = ^(-z); tP± e (z) = ^(z). (15) 

Our aim in this paper is to describe the general eigenfunction of energy E in this 
representation: up to a global (removable) constant factor, it reads as ip E ( z ) = 
ip+(z) + \if) E (z), for any complex projective parameter A, i.e. A G C = CU{oo}. We 
can immediately restrict attention to E > due to the second of eqs. (|l5|). We will 
be particularly interested, on the one hand, in the global profile of these functions, 
and on the other hand, in the position of their zeros because these form the main 
skeleton of the Hadamard product representation (||) (|i~4f . 

The motivation is to better describe the eigenfunctions of a general 1-d Hamil- 
tonian for eigenvalues close to a classical saddle-point energy value. Such an eigen- 
function cannot be simply of WKB form near the saddle-point; instead, it should 
be microlocally modeled by an eigenvector ip E of the dilation operator H near 
(ijP) = with E ~ (up to straightforward displacements, in phase space and in 
energy). Thus, in a Bargmann representation, the eigenfunction near z = ought 
to behave like ^ E {z) for some A, which is the one quantity whose actual value is 
determined by global features of the solution. (In particular, for a parity-symmetric 
1-d Hamiltonian and for a saddle-point located at the symmetry center, only even 
or odd solutions ever come into play; hence parity conservation preselects the two 
eigenfunctions of H with the special values A = +1 and —1 respectively.) 
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3.2 Main analytical results 

The functions ip E (z) are also solutions of the Schrodinger equation written in the 
Bargmann representation, 

l -(-h 2 ^ + z 2 )^ E (z)=E^ E (z). (16) 

It is convenient to use the rotated variables Q, P and Z = 2~ 1//2 (Q — iP) = z e -l7r / 4 
in parallel with q,p, and z. In those variables the quantum Hamiltonian reads as 
the quadratic-barrier Schrodinger operator H = ^(—h 2 d 2 /dQ 2 — Q 2 ). Its Bargmann 
transform happens to be exactly the same operator in the Z variable, simply con- 
tinued over the whole complex plane, so that the eigenfunction equation can also 
be written as 

\{- h2 4pi - Z 2 )* E (Z)=E* E (Z). (17) 



dZ 2 

At the same time, the Bargmann representations obtained from q and Q are equiv- 
alent under a simple complex rotation, 

9(Z) = ij)(z) with Z = ze"™/ 4 . (18) 

Consequently, as a main first result, the above solutions are directly related to 
the parabolic cylinder functions D u (y), defined for example in 17]. As a matter of 
fact, we have : 

•i (Z) = ^ W =*|Jr(i + f ) d_j_*(WI*V (i 9 ) 



Up to rescaling, the Bargmann eigenfunction 4>\( z ) 1S then simply a linear combi- 
nation of known functions, [XD v (y) + D v (—y)] with the notations y = \/2/%z and 



u 2 L h- 



The situation simplifies even further on the separatrix E = 0, where parabolic 
cylinder functions reduce to Bessel functions, as 

D- 1/2 (y) = ((27r)- 1 y) 1/2 ^ 1/4 (y 2 /4) 
D_ 1/2 (y) + D_ 1/2 (-y) = e+-/ 8 (vry) 1 / 2 J_ 1/4 (iy 2 /i) (20) 
-D_ 1/2 (y) + D_ 1/2 (-y) = e^^^y) 1 / 2 J l/A (iy 2 /A) . 

By virtue of eq. (fLG|), eqs. (19) and (|fj) yield the Husimi densities in closed 
form, for all eigenfunctions; e.g., for tp E (z) and ^±i(Z) respectively, 

"?<^> - M^m 10 -^-^ 2 ^ (21> 

H° ±1 (Z,Z) = \K\ 2 J2^/h\Z\ \J Tl/4 {Z 2 /2h)\ 2 e- z ^/ h . 
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This is an interesting extension of earlier similar formulas for the Wigner functions; 
e.g., for ^f, || 



_J_|__ %) lFl (1 + 1; _ 4w/ ,) (22) 

l-K^Trtfa) J (4ipq/h) 



where v\\\x) cx L u (x) (Laguerre functions). These expressions are specially 

simple (constant along connected orbits) because the Wigner representation is ex- 
actly //-invariant under the symplectic evolution generated by the bilinear Hamilto- 
nian H. Inversely, coherent-state representations can never preserve such a dynam- 
ical invariance of hyperbolic type, hence closed-form results like ( |2l| ) are necessarily 
more intricate than their Wigner counterparts and not simply transferable there- 
from. 

Fig. 1 shows contour plots for some of those Husimi densities, all with the 
normalizations K = K = 1. The first example (Tl^(z,z) for E = +1) is plotted 
twice: once with equally spaced contour levels starting from zero (linear scale), and 
once with contour levels in a geometric progression decreasing from the maximum 
(logarithmic scale). The linear plot emphasizes the high-density modulations which 
control the measure concentration of the Husimi density; the logarithmic plot re- 
veals the subdominant structures and especially the locations of the zeros z m which 
determine the Hadamard parametrization. For the parabolic cylinder function itself 
as an example, the factorization formula (||) specifically gives 



Dy{z) 

p(y) 



e p{z) J] ( I 



-) exp (£+^(£) : 



= log D„(0) + (log D u )'(0) y + (log D v )"(0) y 2 /2 



log 



r(^) 



l-u\ 



V2 



7T 



2> 



V 1 

2 + 4 + 



r(-|) 



y 2 . 



(23) 



In order to save figure space, we do not provide the log-plots used to locate the 
zeros for the other Husimi densities but only their linear contour plots, with the 
zeros superimposed as small dots. (The same uniform contour level spacing is used 
throughout to make comparisons easier.) 



4 Asymptotic expansions 

We can then rely upon the known asymptotic properties of the D v {y), which follow 



from the integral representation (14), to investigate two asymptotic regimes for the 
eigenfunctions of H. Firstly, when E/U is kept finite, the asymptotic expansions 
will be valid in the limit y — > oo; this corresponds to eigenenergies very close to 
the classical separatrix energy E = 0. Secondly, if we fix the energy E at a non- 
vanishing value and let fi — > 0, we have to use a different type of asymptotics, 
namely usual WKB expansions. 
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4.1 Energies close to zero 



We use the expansions for D v (y) when \y\ — ► oo [17| at fixed u, which are obtained 
from integral representations using Watson's lemma, and take different forms in 
various angular sectors, 

-^<arg( y )<+| : D v (y) ~ ^e^ £ (Z^ ~ ^/'(l + Q^)) (24) 



n>0 



-7T < arg(y) < -£ : ~ e^l + 0(y~ 2 )) - r j,— 1 e^/ 4 (1 + 0(y- 2 )) 

2 F{-u) 



+ | < arg(y) < +vr : ~ y» e^^l + 0(y- 2 )) - e™ y— 1 ^ (1 + 0{y- 2 )). 

The sectors are specified here as non-overlapping and bounded by Stokes lines, 
i.e. curves of maximal dominance of one exponential factor over the other. (Each 
asymptotic expansion actually persists in a larger sector overlapping with its neigh- 
bours, but this extension will not be of use here.) The above expansions are 
valid for \y\ — > oo within each sector and provide approximations to the shape 
of the eigenfunctions and the positions of their zeros for large |y|; we will use 
them to leading order only (up to 0(y~ 2 ) terms). For the general eigenfunction 
ip^ oc [\D v (y) + D v (—y)\, eqs. (p4|) straightforwardly generate four different ex- 
pansions in the four z-plane quadrants Sj, j = 0,1,2,3 (named anticlockwise from 
S = {0 < avgz < vr/2}). 



4.2 WKB expansions for a fixed non- vanishing energy 

The previous expansions are inapplicable when h — ► with the classical energy kept 
fixed at a non-zero value (in the following we will suppose E > 0). In this regime, we 
have to use WKB-type expansions instead. These can be obtained from the integral 
representations of the solutions fll4|) , by performing saddle-point approximations ; 
equivalently, they can be found directly from the Schrodinger equation written in 
Bargmann variables (|l7|). We will use the Z variable for convenience. 
The general WKB solution can be written, to first order in h, as 

V E (Z) ~ (2E + Z 2 )- l / A (a{%) Q +i<t>{Zo,z)/h + ^ e -i<j>{z ,z)/h^ ^ ^ 

where the exponents are now the classical action integrals, taken from an (ad- 
justable) origin Zq, 

<t>(Z ,Z)= [ Z P E (Z')dZ', P E (Z) = V2E + Z 2 (26) 
JZo 

with the determination of the square root Pe{Z) fixed by the cuts indicated on fig. 
2 (left) and by -Pe(O) > 0. This approximation is valid for h — ► 0, when Z stays far 
enough from the two turning points Z± = ±iy/2E in the sense that \4>(Z±, Z)\ ~^>U. 
The coefficients a and (3 a priori depend upon fi and the region of the complex 
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plane where Z lies. More precisely, the complex Z-plane is to be partitioned by the 
Stokes lines, specified for each turning point by the condition 



i<j>(Z±,Z)/h real. (27) 

Three such lines emanate from every turning point. When the variable Z crosses 
a Stokes line, the coefficients a and (3 change according to connection rules (see 
1 19] for instance); the application of these rules yields the global structure of the 
solution. For full consistency, no Stokes line should link two turning points; this 
restriction forces us to slightly rotate fi into the complex plane, as Ti — * e lt fi, with 
the resulting partition of the plane drawn on fig. 2 (left). The explicit form of 
exponentially small WKB contributions is generally sensitive to the choice e = ±0. 
However, such is not the case for the subsequent results at the order to which they 
will be expressed, so that we may ultimately reset e = 0. 

Before studying a particular solution, we introduce the hyperbolic angle variable 
= arcsinh(Z/ V 2E) which allows to integrate the action in closed form, as 

<f>(Z , Z) = f [sinh(20') + 26'} * = \ \z'^2E + Z' 2 + log(Z' + V2E + Z' 2 )l Z ; 

Z Z L J Zq 

(28) 

the turning points correspond to 9± = ±z7r/2; the action values 0(0, Z±) = ±inE/2 
are frequently needed. 

We fully describe one eigenfunction as an example, \£+(Z) (corresponding to 
A = 0). We identify its WKB form first in the regions 5', <Si, 52, by noticing that 
this eigenfunction must be exponentially decreasing for Z — ► oo in a sector around 
argZ = Sir/A (i.e., z — > — oo) overlapping with those three regions, and then in 
the remaining regions by using the connection rules. The result is 

~ J^^ e+ ^ Z)/H in the regions 5', 5 l5 5 2 

yE( Z) e -,E/2nf e+ i4>{z + ,z)/n +ie -i4>{z + ,z)/n\ [q Sq (29) 

Pe\Z) ' ^ ' 

sSff{Z) ~ ? n) i/ 2 e+7rE/M (e-K^-'^^-i e-W-W) in s 3 . 
Pe{Z) 

The overall normalization factor C(K) is determined by comparison with the direct 
saddle-point evaluation of the integral (|l4|): 

C{%) =K(2irh) l ^e- nE / 4n e- i7T / s (E/e) lE / 2h . (30) 

Eqs. (p^) readily yield the WKB expansions for the general solution as 



Pe(Z) 1 / 2 
C(h) 



n{z) ~ PeJzW 6 +( A - c -) e 



3 +i0(o,z)/^ e -i<t>(a,z)/h 



-i<t>(o,z)/h \ p -i0(o,z)/ft' 



in 5' 



in 5i 



o 
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*f(Z) 
with the notations 



C(h) 



Pe(Z) 1 / 2 

cm 

c{h) 



(1 - Ac+) e +i W> z y n +A e-Morf/n 
e +i<t>(o,z)/n _ x e -mo,z)/h 



■ e ±ivu = ±ie 



±-kE/H 



(c- = l/c+). 



in Si (31) 
in S2 
in S 3 , 

(32) 



5 Large values of the Husimi density 
5.1 In the WKB framework 

We study the particular solution ty+(Z) for a fixed positive energy E as an example. 
From the WKB expansions (29), we derive the Husimi density of this solution, using 
the hyperbolic angle as variable. In the regions 5', S\, S2, away from the turning 
points, we obtain 

\C(h)\* 



n%{z,z) 



[£(cosh 25R(0) + cos29(0))] 1 /2 



exp 



{!(- 



cosh 25ft(0) [sin 29(0) + 1] + cos 29(0) 

(33) 



29(0) 



This formula shows that the Husimi measure concentrates semi-classically along 
the maxima of the exponential factor. Since the variable is restricted to the strip 
|9(0)| < vr/2, those maxima occur on the line 9(0) = — tt/4, which corresponds 
exactly to the branch of hyperbola of energy E in the half-plane 9(Z) < 0. In 
the region Sq D {9(Z) < 0}, H+(Z, Z) also obeys eq. ( |33"|) up to exponentially 
small terms, so that the discussion concerns the whole E'-hyperbola branch in the 
lower Z-half-plane. The above expression simplifies around this maximum curve, 
according to the following remarks. First of all, the variables (Z, Z) are (up to a 
factor —i) symplectic transforms of the original variables (q,p), so the expression 
of the classical energy E = —\Z + V(Z), where V(Z) is analytic, implies the 
following classical velocity along the ^-energy curve : 



dE 



-= = -iZ 



(34) 



and along this curve, we also have \Z\ 2 
pose a small variation 59 as 

5Z\\ = 
5Z\ = 



^cosh 25R(0). Furthermore, if we decom- 



/ 2^cosh0 $1(69) 
v / 2^cosh0 9(<50) 



(35) 



(5Z± is a variation of Z perpendicularly to the E- hyperbola) , then we obtain the 
following expression of the Husimi density around this maximum curve: 



(36) 
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that is, the density decreases as a Gaussian of constant width normally to the 
maximum curve, its height being given by the inverse of the phase-space velocity. 
This corresponds semi-classically to a conserved probability flux along this curve, 
and confirms earlier predictions |HJ (cf. fig. 1, top left and bottom right). 

As a new feature, by contrast, in the upper Z-half-plane there is a maximum 
curve in the region So only, and well below the anti-Stokes line where the zeros of 
ty+(Z) lie. This maximum curve is given by ^s{6) = +vr/4, i.e., it is the other branch 
of the .E-hyperbola. Around it, the Husimi density behaviour is precisely eq. ( |36|) 
times the constant factor exp(— 2nE/h), an exponentially small contribution com- 
pared to that from the lower half-plane; hence this enhancement is semi-classically 
"invisible" but can be guessed on the log-plot in fig. 1, top right. (The correspon- 
dence between the z and Z variables is recalled on Fig. 2, right.) 



5.2 Energies close to zero 

Now using the expansions (^4|) , we can analyze the large values of the Husimi density 
in the case where E/h stays bounded, still in the semi-classical limit Ti — > 0. If we 
still consider the function ip+(z), we find a concentration along an invariant subset 
of the separatrix, i.e., the positive real and the imaginary z-axes. More precisely, if 
\z\ 2 3> h, we have 

Hf{z,z) ~ \K\ 2 V2^K-^- e-2£arg(*)/fi e -29(z) 2 /7i when | arg(z) | < ^ (37) 



\K\ 2 

Hf(z,z) 



cosh(irE/h) 



1 e 2W-,)/R e -2« W Vft when |arg(-*)| < ^. 



We notice that both the longitudinal dependence of the density, and its Gaussian 
decrease away from the separatrix, are exactly the same as for a regular energy curve. 
Thus, away from the unstable fixed point z = 0, the singular limit E — > behaves 
straightforwardly. Here, however, the Husimi density is also described exactly all 
the way down to the saddle-point (for z ~ when E m 0), by the explicit formulae 
@), e.g., 

w ^ o > = Si r (j +! I)f- < 3s > 



We can likewise obtain the rough shape of the Husimi density for a general 
solution As before, whatever type of expansion we use, we find along each of 
the four half-axes (which are asymptotes to the classical energy curves) 

Hf(z,z)~I^e- 2 \ sz ^l n (39) 

where the constant J, depending on A and on the half-axis we consider, can be 
interpreted semi-classically as the invariant intensity of a flux of particles moving 
with velocity z = z along this branch of classical curve. For E 3> 0, the flux 
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is separately conserved along each of the two hyperbola branches: with obvious 
notations, 



1+ = J_j = \K\ 2 V2^h (40) 
I- = I +i = \K\ 2 \\\ 2 V2^h. 

When the energy approaches its critical value 0, the intensities become 



i + = \k\ 2 V^n (41) 
i_ = iK^V^nix] 2 

T+i = _ gJ! - ||A| 2 e ^/fi +2 ^( A) + e~^/ fi } 



2cosh(7r£/fi,) 
IKI 2 



|| A |2 e -7rE/a _ 2 Q(A) +6^}. 



2cosh(7r£/fr) 

Now the flux is only conserved globally: it is easy to check that J+j + I_j = + 



6 Asymptotic study of the zeros 

The asymptotic geometry of the zeros for a general solution can be deduced from 
the following general principles. When a function f(Z) has an asymptotic expan- 
sion (within a sector) combining two exponential behaviours, the function can only 
vanish when both exponential factors are of the same order of magnitude. Thus, 
the two exponents must have equal real parts: this necessary condition defines the 
anti-Stokes lines of the problem in the complex Z-plane. Zeros of the function can 
then only develop along anti-Stokes lines (always in the large- \Z\ approximation), 
and provided both exponentials are present in the given sectorial expansion, with 
fixed non-zero prefactors. 



6.1 Energies close to zero - general case 



Here the expansions (^4|) have to be used, and it is simpler to work in the rotated 
Z variable. Then the exponential factors read as e +lZ I 2 * 1 and e~ tZ / 2h , and the 
anti-Stokes lines on which they balance each other, namely the bisecting lines Lj of 
the quadrants Sj, are simply the real and imaginary Z-axes (fig. 2, right). Now, for 

each fixed j, two independent solutions <&2^ °f e 1- (0) can ^ e specified by imposing 
single-exponential asymptotic behaviours along Lj, as <I>;jP ~ f±\%) e± lZ2 / 2n , A 
general solution ^(Z) is then proportional to A$+' + <&_ for some A £ C, and 
satisfies 

¥(Z) oc Af ( f\z)e +iZ ^ 2n +f^(Z)e- iZ ^ 2n , Z^ooinSj, 
whence the condition ^(Z) = in the sector Sj asymptotically reads as 



(42) 



Z 2 /h ~ (-l) J 2vrm+ilog \kr®(Z) 



for m — > +oo, 



(Z) = -f ( f\z)/f^'(Z). 

(43) 



U), 
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This yields an asymptotic sequence of zeros, as: Z$ ~ e^/s ^2Trmh to leading 
order; thereupon, to the order O(l) included, 

Z^ 2 /fl ~ (-l) J '27rm+«log r (j) (e ijlT / 2 V 2-Kmfi)+i log A+0(m -1 log m), m -> +00. 

(44) 

The final square-root extraction is straightforward, 

ZgWS ~ e«2 + ■° g r,J ' (eiJ ' /2 /g!^ ) + l0SA + O f 

\ 2V27rm \ m-V^ J J 

(45) 

so that the simpler form ( |44| ) will be preferred for further displays of results. 

When the asymptotic analysis concerns a fixed linear combination given as 
[XD u {y) + D v {— y)], A turns into a sector-dependent function A^(A), which is the 
linear fractional transformation induced by the change of basis {D u (y), D u (—y)} — ► 
{<3?^(Z), $_ (Z)}. The equations for zeros like (|44] ) become singular for the two 
values of A which map to = or 00, simply because they yield the pure <£_ or 
solutions which have no zeros (at least asymptotically) in this sector Sj. 

We now list more explicit results. In the sector So, corresponding to {0 < argy < 
+^/2}, 



XD u (y) + D u (e-™ y) ~ (A + e^y" e"^/ 4 1 e+^ 4 . (46) 

Upon the substitutions A<°) (A)r<°) (Z) = (\/27r/r(-i/))y- 2 '' _1 /(* + e~*™), y = 
\J1jTi e* 71 "/ 4 Z and z/ = — \ — ^ , the asymptotic equation ( 44 ) for zeros becomes 

7(0)2 „ + 
"' ~ (2m-l)ir-— log Amiri-i log ^^=M-ilog(\+ie- nE/h ), <argZ<+- 



h v ' h ° V2tt y ' 4 ° 4 

(47) 

A geometrical interpretation will prove useful. Let C± be the two circles in the 
A-plane respectively specified by the parameters (cf. eq. p2|) and fig. 3) 



centers: c± = ±i e ±nE / n , radii: R± = (1 + e ±27r£ / a ) 1/2 (48) 

(they are both centered on the imaginary axis and intersect orthogonally at A = +1 
and —1). Let us also write, for real t, T(^+it) / \/2tt in polar form as (2cosh7rt) _1//2 e J ®(*\ 
defining the phase O(i) by continuity from 0(0) = 0. The formula (|47|) , and its part- 
ners in the other sectors, then become: 

in So = {-vr/4 < argZ < +vr/4}, 

Z (0) 2 E A - c_ 

(2m- 1)tt - — log 4mvr + 9( J B/?i) -« log— (A^{c_,oo}) (49) 

in Si = {+7T/4 < argZ < +3tt/4}, 

7 (1)2 F A- 4 -r 

_ (2m - 1)tt+ - log 4mvr - 0(£?/n) + t log ± (A ^ {0, c_ }) (50) 

n n R+ 
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in S 2 = {-5vr/4 < argZ < -3vr/4}, 

v(2)2 „ A _1 -c 

(2m - 1)tt - -r log 4mvr + 0(£/ft) - % log (A ^ {0, c+}) (51) 

in S 3 = {-3vr/4 < argZ < -tt/4}, 

Z$ 2 ,„ ^ E, _ . A 



(2m - 1)tt + — log 4m7r - 9 (E/H) + i log o ^ (A {c+, 00}). (52) 

U H, ri_|_ 

For general values of the parameters we cannot get more precise information this 
way, except by going to higher orders (but still in the asymptotic sense). We have 
yet no idea about the position, or even existence, of small zeros. Our subsequent 
strategy will be to start from very special cases for which the pattern of zeros is 
well known, and from there to vary continuously the parameters E and A and keep 
track of the zeros along these deformations: zeros are topological defects, so they 
move continuously w.r. to both parameters. We will then exploit symmetries of eq. 
(|T7|), especially the reality of its solutions (a real function of a complex variable t is 



one satisfying f(t) = f(t)). Eq. ( |T7| ) is a differential equation with real coefficients, 
hence it admits real solutions, whose zeros are symmetrical w. r. to the real Z-axis. 
During a deformation of a real solution, zeros could then only leave (or enter) the 
real axis in conjugate pairs, but at the same time the functions considered here 
(solutions of a second-order equation) cannot develop double zeros; consequently, 
each zero is bound to stay permanently real (or nonreal) in the course of a real 
deformation. 

6.2 Even— odd solutions of zero energy 



We use the E = expressions eq. (|20|) in terms of the Bessel functions J ±1 / 4 (Z 2 / 2Ti) 
to view the pattern of zeros more precisely in this particular case. 

We know that, for real [/,, t _/i J^(t) is a real even function having only real zeros 
±j/i,m, m = 1, 2, • • •, with > and j^ m ~ ir(m + f - 3) for large m. This 
translates into the Z variable as follows. 

For A = +1: the even solution (fig. 1, middle right) 

9° +1 (Z) = K(2^/hf' A Z^ J_ 1/4 (53) 

(a real function for K real), is not just even but also invariant under Z — > iZ; all its 
zeros are purely real or imaginary, the m-th positive zero admits the approximation 



Z<®\x=+i ~ \A(2m7r - 371-/4), (54) 



and all other zeros follow by the rotational symmetry of order 4. 
For A = — 1: the odd solution (fig. 1, bottom left) 

V4 1/9 T (Z r 



m\(Z) = ' J 1/4 ^ (55) 
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has exactly the same symmetries for its zeros as the even one (because the auxiliary 
function e -l7r / 4 l (Z)/Z has all the symmetries of ^^(Z)); besides an obvious 
zero at the origin, the m-th positive zero of ^^(Z) lies approximately at 

zW| A= _x ~ ^(2m7r-7r/4). (56) 



6.3 Even— odd solutions of non-zero energy 

If we "switch on" the energy, keeping A = +1, the solutions will not be Bessel 
functions any longer, but they still exhibit interesting features. Equation ( |T7| ) has 
real coefficients, hence its even and odd solutions are real up to a constant factor, i.e., 
up to an adjustment of arg K. For the even solution ^^(Z), the reality condition 
is 



*+i(0) = ( 2h ) m r [l + reaL ( 57 ) 

This even solution can be seen as a real deformation of the solution (j53|). It maintains 
the symmetry w.r. to the origin and to the two Z coordinate axes; only the tt/2 
rotation symmetry is lost for E ^ 0. Due to the two mirror symmetries, as explained 
above, the only possible motion of the zeros during this deformation is a "creeping 
without crossing" along the four half-axes, symmetrically w.r. to the origin. This 
can be checked on the O(l) terms of the expansions (|49|-^), and on the sequence of 
plots: fig. 1 (middle right), fig. 6 (right), fig. 1 (bottom right). At the same time, 
this deformation allows to properly count the zeros at all energies, by continuity 
from E = where eq. (|54|) does count the zeros: in each of the expansions (f49|-|5"2"|) , 



Zm remains the actual m-th zero on the half-axis Lj if the corresponding complex 
log functions are defined at A = +1 as 

log(l-c+) = log|l-c + |-i arctane +7rB//?l , log(l-c_) = log |l-c_|+i arctan e~ 7rE / h 

(58) 

where the arctan function has the usual range (—tt/2, tt/2). 

The same analysis can be performed for the odd real solution ^^(Z) of Ql7D, a 



deformation of the solution (55), for which the reality condition is 



d^ E 
~dZ 



1(0) = K^pj ' (2^)4+1 r (jj + ^ e-/ 4 real. (59) 



6.4 Real solutions 

We now consider more general families of real solutions, i.e., ^ff(Z) = ^f(Z), 6 
C. These exist only for certain values of A, for which we have to adjust K. Since 
(|T^ ) is a second-order equation with real coefficients, it has the real solutions 

qf E (Z) = KOsrf^Z) +m E l (Z)^j , for kEI*, t£RU{oo}=S. (60) 

Under the change of basis {*f X (Z), ^ E i(Z)} — > {D u (y), D u (-y)}, the set {t G 1} 
is mapped to a circle in the projective A-plane, passing through A|t=o = +1, A|( =00 = 
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— 1. The full circle can be determined long-hand using eqs. (57-59), but more easily 
by asking the expansions (|49|), ( |5ll) to yield asymptotically real zeros as reality 
demands: the resulting A-circle is C_ (note that A G C_ <^=^ 1/A G C_). 

Let us now analyze the motion of the zeros in the four sectors as we vary A 
around C_ anti-clockwise from 1 to e 2ln . We already know that the real zeros 
cannot leave the real axis by symmetry. In the expansions ( |49D and ( |5lD for So an d 
S*2 respectively, the only modifications are that arg(A — c_) increases by 2ir, and 
arg(A _1 — c_) decreases by 2ir, inducing a re-labeling of the large zeros in those 
two sectors. Hence each large positive or negative zero creeps to the right until it 
reaches the former position of its right neighbour after one cycle. The small zeros, 
trapped on a real bounded interval in-between by reality, and unable to cross one 
another, can then only follow the same homotopic pattern of behaviour. 

In each of the two remaining sectors, by contrast, the zeros are not confined to 
the imaginary axis whereas the cycle C_ is homotopically trivial in the Riemann 
surface of the relevant logarithmic function. According to the expansions in S± 
(resp. S3), the large zeros in these sectors perform a clockwise (resp. anti-clockwise) 
cycle beginning and ending at their location on the imaginary axis for A = 1. The 
geometric relation A — c+ = R\(X^ 1 — c + )~ l is the asymptotic remnant of reality: 
{fy{Z) = =^ ty{Z) = 0}. This motion in S\ and S3 has been thus shown only for 
large zeros, but we may argue a similar behaviour for the smaller ones by homotopy; 
moreover, by reality these zeros cannot cross the real axis and are confined to a 
bounded region around the origin (contrary to the real zeros, they cannot migrate 
to infinity when A keeps revolving around C_); orbits of non-real zeros during the 
described A-cycle are also symmetrical w.r. to the imaginary axis, thanks to a 
A <-> 1/A symmetry of the real solutions. 

The evolution of all the zeros under the A-cycle C_ is globally depicted for E = 
on fig. 4 (left). 

Similarly, the solutions corresponding to A G C+ can be chosen real w.r. to the 



variable iZ (observing that eq. (17) is invariant under the change {Z — > iZ,E — * 
—E}), and analyzed likewise. 

6.5 Singular values of A 

Each of the expansions (49-152) becomes singular for two values of A from the set 
{0, 00, c_ , c + }. The solutions are then vp±(z),ifj± E (iz), corresponding to "pure D v 



functions", and eqs. ( |24T ) clearly show that D v (y) has large zeros only in two out of 
four sectors, namely along the pair of adjacent anti-Stokes lines {argy = ±37r/4}. 
That is why the zeros in the two other sectors must "escape to infinity" as A moves 
towards one of its special values. Let for instance A decrease from +1 (the even 
solution case) to the special value along the interval [0, 1]. For simplicity we first 
restrict ourselves to the case E = 0, where the eigenfunctions are combinations of 
Bessel functions. The zeros' expansions 

(H|-|5D (with ^±1^=0 = \/2) then become 



singular in the sectors S± and S2, whereas they stay perfectly uniform in the two 
other sectors. In S2, eq. (|5l"l ) holds uniformly as Z^ 2 jfi — > 00, which amounts 
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to — log A + 27rm 3> 1, and the same conclusion is reached for eq. ( pQ) in Si. As 
A — > +0, the two formulae merge into a single one, 

Zl/h ~ i log + (2n - l)vr in Si U S2 for any n G Z (61) 

where the global counting index n S Z matches with m in the sector 5*2 for n S> 
and with (— m + 1) in Si for n <C 0. These zeros thus tend to follow the hyperbola 
branch {2^(Z)9(Z) = ft log X/y/2, R(Z) < 0}, which itself recedes to infinity as 
A — ► +0, as can be seen on fig. 4 (right) followed by fig. 1 (middle left). 

This description can be generalized to the case of a non-vanishing (but small) 
energy. The asymptotic condition for zeros in Si U S2 when A —* must now be 
drawn from eq. ( [43|) itself and reads as 

Z ~ (2n - 1)tt — * log ( A J"^ 1 - f log ( 2^ ) + @(E/h) for n -► ±00. 



(62) 

For small values of A, these zeros tend to follow asymptotically two half-branches of 
hyperbolae which differ as n — > +00 or —00 (because of the Q(logZ^) contribution); 
the matching between this set of zeros for small A and the set of zeros of on 
R_ U iR+ can be done as in the case E = 0. 



6.6 WKB expansions of zeros at a fixed non-vanishing 
energy 

We first consider the particular solution ^^_{Z). In the semi-classical limit of eqs. 
p9|), can vanish only within the regions 6*0 and 6*3 and along anti-Stokes lines, 

now defined by the conditions: 4>(Z±, Z)/h real, respectively. From the complete 
set of anti-Stokes lines (shown on fig. 5 top left, with arg?i = henceforth), the 
presently relevant ones are: in 6*3, the imaginary half-axis below and in 6*0, the 
anti-Stokes line from Z + asymptotic to the positive real axis (fig. 5, bottom right). 
The zeros themselves are given asymptotically by the equations 



<p(Z+, Z) = (+m - l/4)7wr, m G N, m > 1 in S 

4>{Z-,Z) = (-m+l/tyM, mGN, m>l in <S 3 . (63) 

Finally, the expansions 

4>{Z + ,Z) ~ +Z 2 /2 when Z -> +00, , Z) ~ +Z 2 /2 when Z -ioo 

(64) 

restore the former large zero behaviours, rrnr for j = 0,3 respec- 

tively. 

We can also see how the zeros in regions 6*1 and 6*2 go to infinity for a more 
general eigenfunction ^ff(Z), when the parameter A decreases from 1 to along 
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[0, 1] as in subsection ( |6.5| ), but now for a fixed non- vanishing energy. We will then 
be able to compare the results in the two frameworks. 

Using the WKB expansions ( |3l| ) in the different regions of the Z-plane, we study 
the equation ^f(Z) = in each of them |2^]. For A = we found zeros along only 
two anti-Stokes lines. For a general value, there will be zeros in all regions (and, 
inasmuch as A varies independently of h, those zeros are not confined near the above 
anti-Stokes lines). Moreover, we know from symmetry properties that for A = ±1, 
the zeros can only lie on the real and imaginary axes. It is also obvious, from the 
different expansions, that the zeros' pattern depends on the ratio A/c_ = Ac + (cf. 
eq- (|D). 

In the range 1 > A> | c_ | , the differences between the formulae ([n]) for regions 
S' , So, 1S2 are irrelevant as far as the position of the zeros is concerned, so that 
the first one suffices to localize the zeros in those three regions by means of the 
equations 

3f<^(0, Z) = —(ih/2) log A (65) 
5i</K0,Z) = Titt/2 modftTr. 

For A = 1, the zeros are exactly real by symmetry. When A decreases towards |c_|, 
the curve (^) gets deformed towards Z+, keeping the real asymptote at 

both ends. At the same time, the zeros in Si U 1S3 stay along the imaginary axis. 
When A ~ |c_|, there are zeros along all anti-Stokes lines from Z + and along the 
imaginary axis in £3. When A <C |c_|, the zeros along the anti-Stokes lines in 1S0U1S3 
stabilize, whereas the ones in the regions S± U £2 lie along the curve (|65|), which 
recedes to infinity when A — > as in the previous subsection, keeping iR + and R_ 
as asymptotes. All these phenomena appear on the sequence of plots in fig. 5. Once 
more, we can recover the previous asymptotic large zeros by expanding the action 
integrals (f)(0, Z) along the four half-axes. 

7 Conclusion 

The above study should help to better understand more complicated 1-d systems 
when the eigenenergy is very near a saddle-point value and standard WKB theory 
fails. 

This can be illustrated upon the eigenstates of a quantum Harper model. The 
classical Harper Hamiltonian is = — cos(2nP) — cos(2ttQ) on the torus phase 
space (Q mod 1,P mod 1). We have quantized it on the Hilbert space of wave- 
functions with periodic boundary conditions, which has a finite dimension N to be 
identified with (27rh)~ 1 , and taken N = 31 for calculations. We have then selected 
the eigenfunction ^f n immediately below the separatrix energy E = 0; it is an even 
state having the quantum number n = 14, and its Husimi density over the torus is 
plotted on fig. 6, left (cf. also Q 1990, fig. 2a). Now, by expanding the Hamilto- 
nian near one saddle-point like (Q = 0, P = —1/2) we recover a quadratic-barrier 
problem. This suggests to compare ^14 with the even dilator eigenfunction ^^(Z) 
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for E/n = 2ttN\E u \/(2tt) 2 and z = /T x / 2 e i7r / 4 (Q - iP)/y/2 » e i7r / 4 10.(Q - iP) 
(see fig. 6 right). We then observe that not only the high densities of both figures 
around the saddle-points fit very nicely (a result to be expected from semiclassical 
comparison arguments [i"3| , |lO||), but also, more surprisingly, the sequences of zeros 
for the Harper eigenfunction match the comparison zeros as well, and not just near 
the saddle-points but practically all the way out to the extremal points; at these 
points the lines supporting the zeros intersect and the correspondence must fail, but 
this happens much beyond its reasonable range of validity anyway. 
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Figure captions 



Fig. 1. Contour plots of the Husimi densities TLf(z,z) (h = K = 1) for selected 
eigenfunctions ifi^ of the quantum dilator H, using uniform level spacing 0.2 (linear 
plots), except for the top-right plot (drawn in a logarithmic scale using the uniform 
level ratio 1/e from the maximum downwards). The zeros of the Husimi densities, 
constituting the stellar representation, are shown as dots (except in top row). Max- 
ima in the E = plots: H% ^° ax » 3.6963 at z » 0.5409 (at middle left); Tiff (even 
state, at middle right) very steeply peaks at z = 0, with maximum value ~ 10.488; 
Ti^° (odd state, at bottom left) vanishes at z = and reaches its maximum value 
~ 2.0054 at the four points i J zo, zq ~ 1.162. 

Fig. 2. Stokes regions for the asymptotic fi — ► expansions of the eigenfunctions 
of H. Left: fixed E regime (eqs. (p9[) , with argft = e > (the Stokes lines for 
e — ► — e are the left-right transposed of these); Z± are the two turning points; cuts 
for the function Pe{Z) are conveniently placed over the two thicker Stokes lines. 



Right: fixed E/fi regime (eqs. (24),(fj9j)— ( 52)); this picture is the scaling limit of 



the preceding one for E — > when e = (except for the place of cuts) . 

Fig. 3. The complex plane of the coefficient A with the circles C_ corresponding 
to Z-real solutions and C+ corresponding to iZ-ie&l solutions; the finite special 
values 0,c_,c+ are marked (the fourth one being A = oo). 

Fig. 4. Left: positions of the zeros of if)^ =0 (z) for four values of A in anticlockwise 
succession on the circle C_ : Ao = +1, Ai = (y/2—l)i, X2 = — 1, and A3 = (— y/2— 
Right: same for the decreasing positive sequence A n = 10~ n , n = 0,3,6; fig. 1 
(middle left) displays the limiting special value A = 0. 

Fig. 5. Top left: the anti-Stokes lines for the semiclassical expansions of the 
eigenfunctions of H in the fixed E regime. Remaining five plots: positions of the 
zeros of 4>\~ +1 { z )i exceptionally with h = 0.2, for a decreasing positive sequence of 
the type: A = +1; |c_| < A < 1; A = |c_|(= e" 7 ^); < A < |c_|; A = 0. The 
thicker curves are the portions of the H = E hyperbola where the Husimi measure 
concentrates semiclassically. 

Fig. 6. Left: 15 th eigenstate of the Harper model quantized in dimension N = 
l/(2irh) = 31 (with periodic boundary conditions); ^u(Z) has to have a double 
zero at (0, 0) to yield the correct count of 14 zeros within the {Hr < E} region 
(and (31 — 14) outside) |Q; this Husimi density contour plot exceptionally uses 
level spacing 0.089 to match the spacing-to-maximum ratio of the comparison plot 
(next). Right: the Husimi density H+i serving as comparison function, plotted as 
in fig. 1 (maximum value ~ 6.752, at the origin). 
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Stokes lines for E = 1 , e = 0.1 
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Anti-Stokes lines for E = 1.000 E = 1.000 X = 1 
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